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Abstract 

Two chains of ultrasmall Josephson junctions, coupled capacitively with each 
other in the two different ways, straight and slanted coupling, are considered. 
As the coupling capacitance increases, regardless of the coupling scheme, the 
transport of particle-hole pairs in the system is found to drive the quantum- 
phase transition at zero temperature, which is a insulator-to-superfluid transi- 
tion of the pairs and belongs to the Berezinskii-Kosterlitz-Thouless universal 
class. The different underlying transport mechanisms for the two coupling 

schemes are reflected in the difference between the transition points. 
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I. INTRODUCTION 



Systems of ultrasmall tunnel junctions composed of metallic or superconducting elec- 
trodes have attracted considerable interest owing to the significant roles of the Coulomb 
interaction in them. First of all, a sufficiently large charging energy leads to the Coulomb 
blockade effect which exhibits single charge (electron or Cooper pair) tunneling0. In or- 
der for this tunneling to occur, however, it should be energetically favorable with respect 
to the electrostatic energy of the system. Otherwise, more complex elementary processes 
that involve several charge-tunneling events become dominant. In particular, a recent the- 
oretical prediction! and an experimental demonstration! have revealed the cotunneling of 
electron-hole pairs in two one-dimensional (ID) metallic tunnel-junction arrays coupled by 
large inter-array capacitances (see Fig. |l]). Such cotunneling of electron- hole pairs results 
in the remarkable effect of current drag: The current fed through either of the chains in- 
duces a secondary current in the other chain. The primary and the secondary currents are 
comparable in magnitude, but opposite in direction. 

A similar current drag effect has also been observed in a slightly different configuration of 
the system, in which each electrode in one array is coupled aslant to two adjacent electrodes 
in the other array (slanted coupling; see Fig. |^)i. Unlike the case mentioned above (straight 
coupling; see Fig. |^), where the low-energy state can be preserved only if the electron 
and the hole tunnel simultaneously (cotunneling), in this case the low-energy state can be 
preserved for sequential tunneling of the electron and the hole. It has been suggested that 
this correlated sequential tunneling might be more likely than the second-order process of 
cotunneling via a quantum-mechanical virtual state. 

More interestingly, when the tunneling junctions are composed of ultrasmall supercon- 
ducting grains, the counter part of the electron-hole pair becomes the pair of an excess 
and a deficit Cooper pair, which will be simply called a particle-hole pair. Furthermore, 
in such ultrasmall Josephson-j unction systems, the competition between the charging en- 
ergy and the Josephson coupling energy is well known to bring about the noble effects of 
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quantum fluctuationsETO. In a previous worka, it was proposed that, combined with these 
quantum-fluctuation effects, the cotunnehng of particle-hole pairs in capacitively coupled ID 
Josephson junction arrays (JJAs) drives the insulator-to-superfluid transition of the pair^ . 

In this paper, the previous worki on two capacitively coupled ID JJAs is extended to 
consider both straight and slanted couplings. The focus will be on the similarities and 
differences between the two coupling schemes. It is shown that on long-time and long-length 
scales, the two cases are indistinguishable; as the coupling capacitance increases, both of 
them exhibits Berezinskii-Kosterlitz-Thouless (BKT)-type insulator-to-superfluid transition, 
whose superfluid phase is uniquely characterized by the condensation of particle-hole pairs 
and by perfect drag of supercurrents along the two chains. The correlated tunneling nature 
of the elementary process in the slanted case is reflected by its lower transition point. 

Capacitively coupled JJAs can presumably be realized in experiment, even by current 
techniques. Recent advances in microfabrication techniques have already made it possible 
to create large arrays of ultrasmall Josephson junctionsS. Furthermore, the experimental 
realization of the capacitively coupled submicron metal-junction arraysi'S illustrates that 
large interarray capacitances can also be fabricated between two JJAs. In passing, capacitive 
coupling should be distinguished from Josephson coupling (allowing interarray Cooper-pair 
tunneling). The quantum- fluctuation effects in the latter case, as in the case of quantum 
Josephson ladders, have been studied in the literature0. 

The paper is organized as follows: In Section the models for the systems with straight 



or slanted couplings and the regions of interest in parameter space are deflned. Section III is 



devoted to the transformation of the models into equivalent 2D systems of classical vortices. 



In Section the conductivity of the system is examined in the vortex representation. The 
results from Sections |IT| and ^ provide the bases of this work on which the phase transition 
and the current-drag effect are discussed in Section |V|. Section |V| constitutes the main part 
of this paper and present a thorough discussion of the role of particle-hole pairs in the 
quantum-phase transition and the transport in the system. Finally, the paper is concluded 
in Section IVll 



II. MODELS 

Each of the two chains {I = 1,2) of Josephson junctions considered here is characterized 
by the Josephson couphng energy Ej and the charging energies EQ=e^ /2Cq and Ei=e^ /2Ci 
associated with the self-capacitance Cq and the junction capacitance Ci, respectively (see 



which the electrostatic energy Ei=(? j^Ci is associated. Two different ways of coupling are 
considered: Each island in one array is coupled one-to-one to one island (straight coupling; 
Fig. 0) or aslant to two islands (slanted coupling; Fig. ^ in the other array. There is 
no Cooper-pair tunneling between the chains. The intra-chain capacitances are assumed 
to be so small [Eq^ E\ ^ Ej) that, without the coupling, both chains would each be in 
the insulating phaseil. It is also assumed that the coupling capacitance is sufficiently large 
compared with the intra-array capacitances: Ci ^ Cq^C\. In that case, the electrostatic 
energy of the particle-hole pair (~ Ei) is much smaller than that of an unpaired charge 
(~ Eo,Ei). For the most part, this work is devoted to the case of identical chains, but 
non-identical chains will also be briefly discussed. 

In the specified region of parameter space, the time scale of the relevant dynamics is 
determined by the coupling capacitance (C/) and the corresponding Josephson plasma fre- 
quency \^8EjEj / h) . To be explicit, for straight coupling, it is convenient to rescale the 
capacitances by 2Cj (i.e., Co/2C/^Co, Ci/2C/— s>Ci), the energies by hujp = y/AEjEj, and 
the times by 1/ujp. For slanted coupling, on the other hand, a more convenient choice is 
a capacitance scale of 4C/ and a frequency scale of hup = y/2EjEj. Although it is not 
essential, throughout this paper, I adopt these normalizations of the physical quantities for 
simplicity of notation. 

Then, the system with straight coupling can be well described by the Hamiltonian (in 
units of hjp) 



Fig. |I]and Fig. 0). The two chains are coupled with each other by the capacitance C/, with 
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V2KY,Y1 cos[0(£;x)-0(£;x')] 



(1) 



<xx'> 



where the couphng constant has been defined by Ej / 8Ej. The number n{i; x) of excess 

Cooper pairs and the phase (f){i; x) of the superconducting order parameter on the grain at 
X in the chain i are quantum-mechanically conjugate variables: [n{i; x) ,(/){£'; x')]=iSxx'Su' ■ 
The Fourier transform of the capacitance matrix in Eq. (0) takes the following form (in units 
of 2Ci): 

\ 



C(g) = C{q) 



v 



1 -1 
-1 1 



(2) 



where C{q) = Co+CiA(g), with A(g) = 2(1— cosg), is the Fourier transform of the sub- 
matrix C{x,x') within either of the chains. For slanted coupling, on the other hand, the 
appropriate Hamiltonian reads as 



^ = E E ^)^"' ^x^' 



e/' x,x' 



cos[0(£;a;)-0(£;x')]. 



(3) 



<xx'> 



In this case, the capacitance matrix (in units of 4C/) is given by 



C(g) = C{q) 



+ 



1 -1 



-1 



V 



— e 



-iq 



(4) 



III. 2D CLASSICAL VORTEX REPRESENTATIONS 

In this section, I transform each of the models in Eq. (|l]) and Eq. (^ into equivalent 
2D classical system of vortices. The resulting system reveals clearly the nature of the phase 
transition that will be discussed in Section 0. 
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A. Straight Coupling 



It is convenient to write the partition function of the system in the imaginary-time path- 
integral representation as 

^ = n E r ^) [-'^] 

e,x,Tn{e;x,T) 

with the Euchdean action 

s = — ^ E ^(^' ^)^"' ^x^' ^) 

f/' X,X',T 

- V2K Y cos V^0(£; a;, r) 

where V^^ and V,- denote the difference operators with respect to x and r, respectively, 
and the (imaginary-)time slice 6t has been chosen to be unity (in units of l/co'p)Ei. The 
highly symmetric form of Eq. with respect to space and (imaginary) time makes it useful 
to introduce the space-time 2-vector notation r = (x, r) and analogous notations for all 
other vector variables. We then apply the Villain approximationlil to rewrite the cosine 
term as summation over an integer field {m^(£;r)}. Further, with the aid of the Poisson 
resummation formulalil and Gaussian integration, we also rewrite the charging energy term 
as a summation over another integer field {m.r{i; r)} to obtain the partition function 

/oo 
c?0(£;r) exp{-5} (7) 

with 

K 



^' ^')'^--' ['^-'^(^; ^) - 27rm,(£; r)] [V,0(f ; r') - 27rm,(f ; r')] 
+ ^ E t^-'^^^' ^) - . (8) 

£;r 

The variables (f){i; r) and m(£; r) can be conveniently replaced by 0^(r) = </>(!; r) ± 0(2; r) 
and m''=(r) = m(l;r) ± m(2;r), respectively. In this way, one decomposes the Euclidean 
action in Eq. (^ into the sum S = + S~ with defined by 
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r 

Here, the new capacitance matrices C^(x,a;') have been defined by C^{q) = C{q) and 
C~{q) = 1 + C{q). Now, one follows the standard procedure^lil to integrate out {0^(r)}. 
Apart from the irrelevant spin wave part, one can finally obtain the 2D system of classical 
vortices, which is also decomposed into two subsystems Hy = Hy + Hy with 

= y/27r^Kj2^^i^) U^{r-r) v^{r'). (10) 

r,r' 

where the interactions between vortices are defined via their Fourier transforms 

. m , (U) 

A(9)+C±(g)AM 

It is instructive to notice that the vortices represent some sorts of correlations between 
the two chains in the system. In Fig. ^ is represented a configuration of two vortices, 
f(l;r) = +1 and f (2; r') = +1 (r ^ r'). This configuration gives a pair of a vortex 
f^(r) = +1 and an antivortex v^{r') = —1, which tend to be bound to each other. At 
the same time, for v^, it gives two vortices (f+(r) = v^{r') = +1) which are inclined 
to repel each other. For the configuration with a vortex on one space-time layer and an 
antivortex on the other, we have opposite tendencies for and v~. As I will show in 
Section 0, it is the vortices v~ that play major roles in the quantum phase transition of the 
system. Furthermore, the vortices v~ will be shown to be manifestations of the particle-hole 
pairs. 

In Eq. (p!OD, U^{0) always diverges and gives rise to the vortex number equality condition 
^j.i'(l;r) = X]r^(2i'^)' equivalently, the vorticity neutrality condition J2r'^~(^) ~ ^ 
for v~ . A similar neutrality condition, J2r'^~^ ~ should be satisfied unless Cq = (see 
Ref. More importantly, it should be noticed that the two fields and f~ actually are 
not independent of each other since m^(l;r) and m^(2;r) in Eq. (^, and hence f(l;r) = 
[f+(r) +v^{r)]/2 and v{2; r) = [f+(r) — v^{r)]/2, can take only integer values. As depicted 



with open circles in Fig. ^, {v^,v~) at each r can take only half of the elements in the 
product set of integers ZxZ; t>+ and v~ are topologically coupled with each other. However, 
this topological coupling is irrelevant and can be safely neglected, which will be discussed 
in more detail in Section |V]. 



B. Slanted Coupling 



In the same manner as that leading to Eq. (^, from the path- integral representation of 
the partition function Eq. (^, one can obtain the partition function for slanted coupling: 



/oo 
rf0(£;r) exp{-5} 
. ., , -oo 



(12) 



';r m(^;r) 



with 



S = K J2 C{i,i']X,x')Srr' [V,(^(£;r) - 27rm,(£;r)] [V,(/>(£';r') -2vrm,(f;r')] 

£/';r,r' 

+ i^^[V.0(£; r)-27rm,(£; r)]^ 



(13) 



Unlike the previous case, owing to the last term in the capacitance matrix C in Eq. (0), the 
replacement of the variables (f){i; r) and m(£; r) by 0^(r) and m^(r), respectively, provides 
no help at this stage. Instead, we directly integrate out 0(£; r) to get the Hamiltonian for 
2D classical vortices: 

Hv = 47T^K J2 ^(^; r)[/o(r, r')f (£; r') 

e;r,r' 

+ Att^K J2 v'ir)Uj{r, r>-(r') 

r,r' 

+ An'Kj2y7iir)Uiir,r')v;,ir'), (14) 

r,r' 

where the vortex interactions are given via the Fourier transforms 



Uiiq) 



1 A(g) 
r(q) 4 



C{q) A{q) + \l + C{q)\A{uj) + -A{q)A{u;) 



(15) 
(16) 



with 



T(q)= A(g) + C(g)AH A(g) + 1 + 
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(17) 



In Eq. (p^, the vortices v~ and are defined by the direct differences = [f (1; r) — t>(2; r)] 
and by the slanted differences vJi = [v{l;r) — 1>(2; r+e^:)] (see Fig. The appearance of 
the additional term in f ^ in Eq. ( p!^ is no surprise since with a vortex on one space-time 
layer the correspondent vortex has two alternative nearest-neighbor sites to sit on on the 
other space-time layer. However, at the long times and lengths that we are interested in, any 
configuration of gives the same energy as the corresponding configuration of v~ , which 
can be seen in Fig. ^. Thus, in the last additional term in Eq. (|I^ can be replaced 
simply by f ~, leading to 



Hv ~ ^n^K J2 ^(^; r)f/o(r, r')v{i; r ) 

l;r,r' 

+ Sn^Kj^^'i^Wiir, r>-(r'). 



(18) 



This point can also be seen in a more rigorous way by rewriting the Hamiltonian in 
Eq. (|l^) in terms of f ^(r); i.e., Hy = Hy + Hy + SHy, where Hy has the anticipated form 



= 2n^K J2 (r) (r , r') (r') 

r,r' 



(19) 



with the vortex interactions 
1 



t/+(q) 



T(q) 
1 



C{q) A{q) + {l + C{q)\A{uj) + -A{q)A{uj) 



1 



The term SHy, defined by 



1 + C{q) \ \ A(g) + C{q)A{uj) \ + -A(g)A(^) 



(20) 
(21) 



5Hy = -K J2 E^"W'^^"''^V) 

a,f3=i: r,r' 



(22) 



with 



A{q) 
T(q) 



+A(g) +2ising 
^ —2i sin q —A{q) 



(23) 



describes the interaction between f + and v~ . Since the numerator in the interaction (5f/"^(q) 
is the third order in q, at long times and lengths (q ^ 0), 5Hv can be ignored compared 
with Hy. Moreover, in the low frequency and momentum limit (neglecting the terms of 
order 0{c^) or higher, the vortex interactions in Eqs. (^) and (^) are simply reduced to 
those in Eq. ([TTD for straight coupling. For this reason, it is concluded that, at long times 
and lengths, the 2D vortex representation for slanted coupling is equal to that for straight 
coupling except for the different coupling constant: K/\/2 for straight coupling and K for 
slanted coupling. 



IV. LINEAR RESPONSE 

The mathematical mapping in the previous section enables us to examine the existence 
and the universal class of the phase transition, and yet we need to identify and characterize 
the phases on both sides of the transition (see Section 0). A common method is to measure 
the response of the system to an external perturbation. In this section, I consider the 
conductivity, specifically, the linear response ai^i^iuj) of the current in chain £i to the voltage 
applied across chain ^2 (see Fig. ^ and Fig. ^ and rewrite it in the vortex representation in 
accordance with the previous section. 

The response function ai^^^iui) can be obtained via the analytic continuation 

(^iii2{^) = —lim.gi^i^{q,iuj'^uJ + iO^), (24) 

where Geie2 is the Fourier transform of the imaginary-time Green's function 

ge,e,{x,T) = (T,[J(^i;x,r)J(£2;0,0)]) (25) 

with the time-ordered product Tr and the current operators x) = sin Vx0(^; x). Owing 
to the symmetry between the two chains, the conductivities an and (T21 can be written as 

crn{uj) = ^Mu) + cr4u)], (26) 
a2iiuj) = ^[a4uj)-a4uj)] (27) 
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in terms of a± defined in a manner analogous to Eqs. and (^) with I^{x) = I{l',x) ± 
I{2;x). In the same manner as in Section y, one can get the vortex representation of the 
corresponding Green's functions Q± (see Appendix^): 



(28) 



where 7 = a/2 for straight couphng or 7 = 1 for slanted coupling, and the average {■ ■ ■)y 
is with respect to the total vortex Hamiltonian Hy = Hy + Hy. In Eq. (^8]), the same 
Hamiltonians Hy as in Eq. (|10D and the same interactions t/^ as in Eq. ([TTl) are used both 
for straight and slanted coupling. This should be valid at the long times and lengths we are 
interested in, as discussed in the previous section. One of the benefits of the representation 



in Eq. is that the vortex contribution of the second term in Eq. 28 can be estimated in 



the standard renormalization group (RG) approach for 2D classical vorticeg3il3. 

V. QUANTUM PHASE TRANSITIONS 

Now, I turn to the main subjects of this work, the quantum-phase transition and the 
current-drag effect in the system based on the 2D classical-vortex representations of the 



partition functions in Eq. ( |10D and in Eq. (|T9|) and the response functions in Eq. In 



Section |T|, it was established that, aside from the differences in the coupling constant, the 
two coupling schemes are equivalentat at long times and lengths. First, the discussion will 
focuse on straight coupling, the results of which can be applied in full to slanted coupling 
with the coupling constant properly replaced. Some important differences between the two 
coupling schemes will be discussed at the end of the section. 

It is not difficult to understand separately the physics described by each of the Hamil- 
tonians Hy in Eq. ([T0|). Unless Cq = 0, the length-dependent anisotropy due to C^{q) = 



C(g) ^ 1 in t/+(q) in Eq. (|TT]) fades out on (space-time) length scales larger than \/Ci/Cq 
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and thereby C/+(r) = 27r[C/+(0) — U'^{r)] is simply reduced to the isotropic logarithmic 
interaction 

C/+(r) ~ ^O^logr, (r > 1). (29) 
This results in the usual vortex Hamiltonian 

//+ ~ -7rX^^t;+(r)log|r-r>+(r) (30) 

r,r' 

with the effective coupling constant K'^^ = ^ Ej/16Eo. Because we assumed at the begin- 
ning of the paper that -Bo, -Bi ^ Ej, Kgff is substantially smaller than the BKT transition 
point Kbkt — 2/7r; the vortices v'^ always form a neutral plasma of free vortices regard- 
less of K (i.e., regardless of Cj). In the case of Cq = 0, U~^{r) is only short ranged: 
C/+(r) ~ \/Clexjp{—r/\/C[). In this case, the vortices v+ even form a non-neutral plasma of 
free vortices. In any case, the system of vortices v'^ always forms a plasma of free vortices, 
regardless of the value of K. On the other hand, neglecting the short-distance anisotropy, 
U~{r) is isotropic in space-time and logarithmic in distance: C/~(q) ~ l/[A(g) -|- A(a;)] and 
U (r) ~ log r leading to the 2D Coulomb gas Hamiltonian for v~ 

Hy ~ -V27r^X^v-(r)log|r-r>-(r'). (31) 

r,r' 

It follows that the system of vortices v~ exhibits a BKT-type phase transition at K = 
such that K-/V2 = Kbkt ^ 2/7r. 

At this point, one might be tempted to conclude that, as K is decreased, the total system 
Hv goes through a BKT-type transition at \/2Kbkt which is entirely driven by the vortices 
v~ , with playing no role. This scenario, however, should be carefully tested against the 
topological coupling between v'^ and v" discussed in the previous section. 

For this goal, it is convenient to consider the subsystem {v^} of {v^} satisfying v^{r) — 
for all r. In this subsystem, v*{r) can take only even numbers, and hence the Hamiltonian 
can be written as 
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H*y ^ -7r(4irj,) J]K(r)/2]log|r-r'|K(r')/2] (32) 

r,r' 

unless Co = 0. Under the assumption Eq, Ei ^ Ej {'^K^g ^ Kbkt), there are a significant 
number of free vortices of v*. Obviously, it is also true in the case of Cq = 0. These free 
vortices of v* substantially affect the topological coupling and eventually make it irrelevant: 
As illustrated in Fig. ^, the vortex- antivortex pair of f ~ is always accompanied by two 
vortices of (enclosed in dotted ellipses). The interaction between these two vortices of 

is completely screened out by the free vortices of v* (indicated by double arrows in the 
figure). The two vortices of v^, therefore, cannot affect the interaction energy of the vortex- 
antivortex pair of v^; they can only slightly change the fugacity of f ~. It is concluded that, 
as K is increased, the system exhibits a BKT-type transition at Kf ~ \/2Kbkt — 2"\/2/7i", 
which is exclusively attributed to the vortices v~ . 

To examine the states of the system on both sides of K^*, it is first noted that the 
vortices v^{r) are topological singularities in the 2D space-time configurations of the phases 
(p'^lr) = 0(1; r) ± 0(2; r). The plasma of free vortices of leads to complete disorder of 
0^. Accordingly to the uncertainly relation between 0"^ and the conjugate variable n'^{x) = 
[n(l; x) +n{2; x)]/2, i.e., A0+Ar?,+ > 1, this means that n~^{x) for each x takes well-defined 
value, which should be zero due to the large single-charge(Cooper pair) excitation energy 
of order of Eq or Ei. n~^{x) = 0. Consequently, the variable n~{x) = [n{l]x) — n{2;x)]/2 
conjugate to 0~ represents the number of pairs of excess and deficit Cooper pairs (i.e., 
particle- hole pairs). Further, it is also noted that the effective model in Eq. (|I^) or Eq. ( pi]) 
is a vortex representation of the quantum phase mo delii 

^ X ^ X 

Therefore, the BKT-type phase transition at driven by is nothing but an insulator- 
to-superfiuid transition of the particle-hole pairs: Although particle-hole pairs are always 
the lowest excitations, below Kf, they cannot move along the system without external bias 
due to the Coulomb blockade associated with the charging energy Ej. For K > Kf, on 
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the other hand, the particle-hole pairs condensate to form a superfluid and can move free 
along the system (see Fig. |^. The formation of bound dipoles of vortices v~ is an effective 
manifestation of the condensation of the particle-hole pairs. 

Such particle-hole transport can be confirmed by examining the current responses the 
two chains given by Eq. (^). Due to the free vortices of , it follows directly from Eq. ( P5[ ) 
that cr+(u;) vanishes {uo <^ 1). On the other hand, the tightly bound vortex-antivortex pairs 
of f ~ above K^^ result in the response function 

a_{uj) = ^^(2-^\5{u) (cu<l), (34) 



V Kr 
and hence 



aii(^) = -a2i(^) = ^^^^ 



where the renormalized coupling constant Kr is defined by 



2 



^-^-yElrr(-"(rK(0)>- (36) 

^ r 

Thus the system exhibits superconductivity and carries currents along the two chains equally 
large in magnitude but opposite in direction. This perfect drag of supercurrents reveals 
that the charges indeed transport in the form of particle-hole pairs, which are bound by 
the electrostatic energy Ej associated with Cj. For K < K^*, on the other hand, the 
system displays insulating particle-hole I-V characteristics, qualitatively the same as those 
in Refs.i^. 

The argument so far also holds for slanted coupling if only one replaces K/ -\/2 by K; 
the system with slanted coupling exhibits a BKT-type transition at K^'' ~ Kbkt, and the 
superfluid state is characterized by the response functions 

aniuj) = -a^iH = ^ U - <5(^), (37) 



where 

Kr- K 2 



^ = ^-tEM'(^~W^"(0)>- (38) 
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It is interesting, however, to notice that Kf is quite a bit smaller than K^^ (see Fig. 
This reflects the difference between the two coupling scheme in the underlying transport 
mechanism: The correlated sequential tunneling of particle-hole pairs, a first-order process, 
is more likely than cotunneling, a second-order process. 

VI. CONCLUSION 

In conclusion, the properties associated with particle-hole pairs in two capacitively cou- 
pled Josephson-j unction chains, considering both the straight and the slanted couplings have 
been investigated. In particular, the transport of particle-hole pairs was found to drive the 
BKT-type insulator-to-superfiuid transition with respect to the coupling capacitance, re- 
gardless of the coupling scheme. The superfiuid phase {K > Kc) is uniquely characterized 
by the absolute drag of supercurrents along the two chains. 
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APPENDIX A: 

In this appendix, the vortex-representation of the imaginary-time Green's function in 
Eq. p8| ) is derived. For simplicity, only the Q- for straight coupling is derived here. The 
application of the same approach to Q+ should be straightforward, though. In addition, at 
long times and lengths, the derivation should also hold for slanted coupling if one replace 
K/^hy K (see Section fT]). 

In the imaginary-time path-integral representation, the Green's function ^_(ri,r2) = 
(/_(ri)/„(r2)) can be written as 
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where the Euchdean action iS is given by Eq. (|^). By changing the variables from (f){i;r) 
and m(£; r) to 0^(r) and m'^(r), respectively, one obtains 

5^ has been given in Eq. (^. The 0^-integration can be performed easier by first introducing 
an auxiliary field as follows: 

^-(ri,r2) = --^n E / / rf'J°(r) J-(ri)J-(r2)exp{-5+-5-}, (A3) 

where 

V J. J./ V J. J./ 

+ i^J(r) ■ [V0^(r) -27rm^(r)] . (A4) 

r 

Now, integrating out 0^(r), one finally gets the vortex- representation of the Green's function 

(A5) 

where the average (■ ■ ■)y is with respect to the total vortex Hamiltonian Hy = Hy + Hy. 
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FIG. 1. Schematic picture and equivalent circuit of the system with straight inter-chain cou- 
pling. After Averin, Korotkov, and Nazarov, Phys. Rev. Lett. 66, 2818 (1991). 
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FIG. 2. Schematic picture and equivalent circuit of the system with slanted inter-chain coupling. 
After Delsing, Haviland, and Davidsson, Czech. J. Phys. 46, 2359 (1996). 
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FIG. 3. Displacement vortex v~. The configuration of one vortex v{l;r)=+l on space-time 
layer 1 and another v{2; r')=-|-l on layer 2 (r 7^ r') corresponds to a pair of displacement vortices, 
a vortex v~ {r)=+l and an antivortex v~ {r')=—l. 
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FIG. 4. Topological coupling of v'^ and v . At each space-time position r, {v'^,v ) can take 
only half of the elements in ZxZ as depicted with the open circles in the figure. 
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FIG. 5. Displacement vortex for slanted coupling. At long (space-time) lengths, the vor- 
tex-antivortex pair of v~i gives the same energy as that of v~. 
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FIG. 6. Due to the topological coupling, the vortex-antivortex pair of the displacement vortices 
v~ (indicated by dotted ellipses) induces two vortices of The interaction between these two 
vortices is, however, completely screened out by the free vortices of with ±2 (indicated by the 
douple arrows). 
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FIG. 7. Phase diagrams of capacitively coupled Josephson-j unction chains (a) with straight 
coupHng and (b) with slanted coupling. The intra-chain coupling constants Kq and Ki are defined 
by Kq = \/Ej/%Eq oc \/Cq and by Ki = y/Ej/SEi oc respectively. In the text, only the 

region of K S> i^o, -f^i, where particle- hole pairs are relevant, is considered. 
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